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ABSTRACT 

The prospects for the measurement of the strong coupHng 
constant ay^{Mz) to a relative uncertainty of 1 % are dis- 
cussed. Particular emphasis is placed on the implications re- 
lating to future High Energy Physics facilities. 

I. INTRODUCTION 

Quantum Chromodynamics (QCD), the theory of the strong 
interaction, has a single free parameter, the strong couphng as- 
The coupling depends on the renormalization scheme and the 
energy scale, Q. Once asiQ) is determined from an experi- 
mentally measured process, any other process mediated by the 
strong interaction can be calculated to arbitrary accuracy, at 
least in principle. Most determinations of as are based on per- 
turbative QCD, where it is conventional to evaluate the coupling 
in the MS scheme, which is only defined in perturbation theory. 
Furthermore, it is also customary to choose the mass, Mz, 
as the reference scale. We shall adhere to these conventions and 
quote, for the most part, ay^{Mz) in our discussions. 

A precise measurement of is motivated by a number of 
considerations: 

1. The couplings of the electroweak theory, a^m and sin^ 9w, 
have been determined with a precision of about 0.1 %. In con- 
trast, the strong coupling is presently known only to about 5 %. 
It is pertinent to improve the accuracy with which the strong 
coupling has been measured in order to place it on a more equal 
basis with respect to the other interactions. For example, the 
current accuracy of as measurements is one of the main limita- 
tions on Standard Model electroweak tests at LEP and SLC [|l]]. 
In addition, attempts to constrain Grand Unified models, from 
the convergence of the standard model SU(3), SU(2) and U(l) 
couplings at a Grand Unification Scale, are similarly limited by 
the accuracy with which as has been measured. 

2. QCD with its one parameter, as, must account for the rich 
phenomenology which is attributed to the strong interaction, in- 
cluding perturbative and nonperturbative phenomena. A funda- 
mental test then of QCD is the determination of as from exper- 
imental measurements which probe complementary processes. 
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This test is only meaningful if the values of 05 being compared 
have been measured with similar, good accuracy. 

3. The QCD /3-function (which is known to three loops in the 
MS scheme) determines the evolution of the coupling. Accurate 
measurements of 0:5 over a wide range of momenta provide an 
additional fundamental test of the theory. Tests of the QCD (3- 
function constrain physics beyond the standard model, in partic- 
ular models with additional colored particles. Measurements of 
the energy dependence of observables in a single experiment, 
such as jet variables at e+e^ or pp colliders, can also test the 
QCD /3-function. 

The last two reasons given above for an accurate measurement 
of as emphasize that it is not sufficient to determine as using 
a single method, but that precise measurements are necessary 
using different processes and widely different values. 

For the presentation here, we consider the prospect to mea- 
sure ay^{Mz) with 1 % accuracy. We attempt to identify those 
methods which offer the greatest potential for such precision. 

Figure 1 presents a summary of the most accurate mea- 
surements of as which are currently available. Measurements 
performed at scales different from A/^ have been evolved 
to Q^=M| using the three loop QCD /3-function (in the MS 
scheme). All determinations of as receive contributions from 
theoretical systematic errors. These are, in many cases, the 
dominant sources of uncertainty. In general, they are diffi- 
cult to estimate. In determinations based on perturbative QCD, 
sources of such errors are the truncation of the perturbative se- 
ries and nonperturbative effects (such as hadronization). Most 
of the perturbative calculations have been carried out to next- 
to-leading order (NLO), and, in a few cases, to next-to-next-to- 
leading order (NNLO). The sources of theoretical uncertainty 
in determinations based on lattice QCD are discussed in sec- 
tion III. In a few cases, as results are limited by experimental 
uncertainties. 

Theoretical uncertainties are in general not gaussian- 
distributed and are estimated from a variety of different meth- 
ods. Consequently, the correlations between different 0:5 mea- 
surements are difficult to estimate. Given this difficulty, there is 
not a unique procedure to define a world average for the results 
shown in Figure 1 . A number of proposals for world averages 
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Figure 1: Summary of current accurate measurements of as, 
by technique. The as measurements are based on perturbative 
QCD, except where otherwise noted. 



exist [|, 111. We shall use [|] 

ay^{Mz) = 0.118 ±0.005 
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as our nominal value for the world average. 

In conducting this study, we considered a wide range of ap- 
proaches to the measurement of aj^{Mz), eventually identi- 
fying four methods which exhibit the potential to yield results 
with about 1 % precision. The existence of perturbative cal- 
culations to at least next-to-next-to-leading order is a prereq- 
uisite for 1 % accuracy. A number of such calculations are in 
progress [^, and we shall assume that they will be available for 
the experimental measurements in question. 

The four methods are: (1) the evolution of the parity vi- 
olating structure function xF^, (2) the Gross-Llewellyn-Smith 
sum rule, (3) spin averaged splittings in the T and systems, 
and (4) hadronic observables in e+e^ annihilations. Items (1) 
and (2) in the above list are measured in deep inelastic neutrino 
scattering experiments. Item (3) is based on lattice QCD, all 
other methods use perturbative QCD. In addition to these four 
approaches, we found that two other methods, the evolution 
of the parity non-violating structure function F2 at high x, and 
the jet Et spectrum in high energy proton-(anti)proton colli- 
sions, offer the possibility to determine as with good accuracy 
in regions of which are complementary to those of the other 
measurements. The ultimate accuracy with which as can be de- 
termined using these last two techniques is uncertain at present, 
however Therefore we do not include them in our final list of 
techniques which might yield an as result with 1 % precision. 

The remainder of this report is devoted to a presentation of the 
various methods we considered for an ^5 measurement, with a 
particular emphasis on the implications for future High Energy 
Physics facilities. 



Measurements of nucleon structure functions from the deep 
inelastic scattering (DIS) of a lepton on a nuclear target have 
yielded some of the most precise results for the strong coupling 
constant as- The field of nucleon structure functions remains 
very active, with major experiments in operation at CERN, 
DESY, Fermilab and SLAC. New structure function data, ex- 
tending the measurements to previously unmeasured regions of 
the kinematic variables x and and utilizing polarized targets 
and probes, have recently become available. These programs 
are expected to continue for at least the rest of this decade. 

In this section, we assess which of these new data have the 
potential to yield an measurement with 1 % accuracy. We do 
not review the formalism of structure functions or provide more 
than an indication of the methods used to determine as from 
them. References to existing literature with such information 
are given below where deemed appropriate. 

A. DIS Nucleon Structure Functions 

The basic kinematic variables of DIS are the of the in- 
teraction, given by the difference in 4-momentum squared be- 
tween the outgoing and incoming leptons, and the Feynman 
variable x defined by x=Q^ / {2M{E-E')), where M is the mass 
of the target nucleon, with E and E' the energies of the initial- 
and final-state leptons, respectively, as measured in the labora- 
tory frame. In the quark-parton model, x is interpreted to be 
the fraction of the nucleon's energy carried by the struck par- 
ton. Experiments in DIS measure the energy and scattering an- 
gle of the final-state lepton and/or recoiling hadronic system. 
The lepton probes are either electrically charged (electron e and 
muon probes) or neutral (neutrino v or antineutrino v probes). 
The dominant mechanism for charged lepton scattering is sin- 
gle photon exchange in the t-channel between the lepton and 
nucleon system, while that for vovT' scattering is single 
exchange. For values which approach or exceed Mw, 
and Z'^ exchange become important for charged lepton scatter- 
ing. Nonperturbative QCD corrections to single-parton scatter- 
ing contribute higher twist terms to the cross sections, which 
scale like (1/Q^)" (n=l,2,3- • •) and are important at low 
(typically Q^<4-5 GeV^). 

Of the many structure functions necessary to describe DIS 
cross sections in their most general form, only three are can- 
didates for a precise measurement of as- the structure func- 
tions F2{x,Q^), F^iXjQ^) and gi(x,Q^). F2 is measured from 
the neutral current cross section for unpolarized charged lep- 
tons to scatter from unpolarized targets and from the sum of the 
charged current cross sections for neutrinos 1/ and antineutrinos 
77 to scatter from unpolarized targets. F3 is measured from the 
difference between the charged current v and 77 cross sections 
for scattering from unpolarized targets, gi is measured from 
the asymmetry in the cross sections for longitudinally polarized 
charged leptons to scatter from polarized targets if the beam and 
target polarizations are parallel, compared to the case that they 
are antiparallel, and from the corresponding asymmetry for tar- 
gets which are polarized perpendicular to the beam polarization 
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directions [^] . The techniques that have been used to determine 
as using F2, F3 and gi are mentioned in the following section. 

Structure functions can be resolved into color singlet and 
color non-singlet components. In QCD, the singlet and non- 
singlet terms evolve differently with The singlet com- 
ponent receives a contribution from gluon splitting into pairs. 
As a consequence, the evolution of the color singlet term 
depends not only on the running coupling constant asiQ^) but 
also on the probability for gluon splitting, given by the gluon 
distribution function g(x,Q^). This dependence on g{x,Q^) is 
not important if x is larger than about 0.25 because the prob- 
ability for gluon splitting at large x is small. Gluon splitting 
does not contribute to the non-singlet component of the struc- 
ture functions: the evolution of this term depends on as(Q^) 
only, irrespective of the x range. Depending on the nature of the 
target (e.g. deuterium d2 or hydrogen h2), F2 is either a pure 
singlet or a mixture of a singlet and a non-singlet, whereas F3 
is always a pure non-singlet. 

B. Methods used to Determine as from DIS 
Structure Functions 

The following methods have been used to determine as using 
the F2, F3 and gi structure functions: 

1. the evolution of F2 at high x values, measured in 
charged lepton scattering (4 %) [|^; 

2. the same method as given in item 1, except at low x values 

(11%) [§; 

3. the evolution of F3 multiplied by the kinematic vari- 
able X, i.e. the evolution of xF^ (5 %) [p^; 

4. the Gross, Llewellyn-Smith (GLS) sum rule, based on F3 
at fixed Q^, integrated over all x values (7 %) JTl]]; 

5. the Bjorken sum rule, based on the difference between the 
gi structure functions of protons and neutrons at fixed Q^, 
integrated over all rvalues (4%) ljl|]; 

6. the shape of F2 from charged lepton scattering in the limit 
of very large and very small x values (9 %) 1 13 1. 



The reference given after each item refers to the most pre- 
cise result available for the method. This precision itself, 
Aay^{Mz)/a^{Mz), is given by the number in parentheses, 
where ayf^{Mz) is the value of as after it has been evolved to 
the scale of the Z° mass and Aay^{Mz) is the corresponding 
uncertainty including statistical and systematic terms. 

In addition to the methods listed above, DIS experiments have 
measured as using one technique which is not based on struc- 
ture functions: the measurement of as using jet rates |jl4|]. This 
method is very similar to the one based on event shapes from 
e+e^ annihilations and has similar sources of systematic un- 
certainty. The present accuracy of the result for a^{Mz) from 
DIS jet rates is about 8%. It is not likely that this method will 
yield a result for a-^{Mz) with precision better than about 5 % 
unless a next-to-next-to leading order QCD calculation becomes 
available. The overall situation for this measurement is similar 



to that for jet rates from e+e~ collisions and we will not discuss 
it further 

From the above list, it is seen that the most precise DIS re- 
sults for as are obtained from the evolution of F2 at high 
X (item 1), the evolution of xFs (item 3), the Bjorken sum 
rule (item 5), and - with somewhat less precision at present - 
the GLS sum rule (item 4). It is of note that all three structure 
functions F2, F3 and gi contribute at least one measurement 
with 4-5 % accuracy, illustrating the strength of the complemen- 
tarity offered by the unpolarized charged lepton, neutrino, and 
polarized charged lepton programs. The results utilizing the 
evolution of F2 at low x (item 2) and the shape of F2 (item 6) 
are less accurate. Method 2 is unlikely to provide a precise re- 
sult for as in the future, since the evolution of the singlet 
component at small x depends on the gluon distribution func- 
tion g{x,Q^), as mentioned above: this situation will not change 
for future data sets. If data are collected using different nuclear 
targets so that the singlet and non-singlet components of F2 can 
be separated, the evolution of the non-singlet component of F2 
at relatively small x values could still be a viable method for 
an accurate as result: this was not found to be the case in [||, 
however, which included such an analysis. It is more difficult 
to assess the future status of the as result based on method 6 
since this method has only recently been proposed. This method 
is based on the asymptotic behavior of the QCD resumed pre- 
diction for F2 at large and small x and has been applied 
to HERA data. The dominant uncertainty arises from the am- 
biguity in the choice of the renormalization and factorization 
scales [pj|]. This suggests that a reduction in the uncertainty of 
ay^{Mz) below the 5 % level will require the inclusion of sub- 
leading terms, the prospects for which are unknown. Further 
theoretical understanding of this method will probably be re- 
quired before it can be used to accurately measure as- We will 
not consider this method further The remaining discussion on 
the prospects for a precise 0:5 measurement from DIS therefore 
concentrates on items 1, 3, 4 and 5 in the above list. 

C. Future Prospects for a Precise as 
Measurement from DIS 

The future facilities which we consider for the purpose of 
evaluating the potential for a 1 % measurement of a-^{AIz) 
from DIS are the following: 

1 . HERA with a luminosity upgrade, able to deliver data sam- 
ples of about 150 pb^^ per year, yielding a total data sam- 
ple for the HERA experiments of 500-1000 pb^i; 

2. an electron-hadron collider utilizing the LHC, referred to 
as "LEPxLHC" in the following; 

3. a (F) beam from the Tevatron with upgraded luminosity, 
i.e. Tevatron "Run 2" and TeV33, available for fixed tar- 
get experiments; it should be emphasized that the prospec- 
tive fixed target neutrino facility under consideration here 
would make use of the full energy Tevatron beam; 

A. a.v (V) beam from the LHC, available for fixed target ex- 
periments; 
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5. future experiments measuring the polarized structure func- 
tion gi . 

The facilities listed above have often been presented as natural 
extensions to the HERA, Tevatron, and LHC programs, with the 
possible exception of a fixed target facility at the LHC. It is not 
clear whether it is feasible to incorporate a fixed target neutrino 
facility into the LHC program. 



We next discuss each of items 1, 3, 4 and 5 from section |I.B 
in the context of these future facilities. 



Evolution of F2 at High x 



A measurement of as from the evolution of F2 is best 
performed using charged lepton scattering on unpolarized tar- 
gets (for a review of this method, see [|l5[|). The relevant fu- 
ture facilities for this measurement are an upgraded HERA and 
LEPxLHC. 

An as measurement based on the evolution of F2 is not possi- 
ble using the current HERA data sample because of the scarcity 
of data at x values above about 0.20. HERA has kinematic ac- 
cess to X values up to about 0.50 for values greater than 
10^ GeV^, however For the purposes of an as measurement, 
HERA data at 10"* GeV^ and x ^ 0.50 are interesting be- 
cause the large x value ensures suppression of the contribution 
from the gluon distribution function, while the value is sim- 
ilar to that in Z*^ decays: this offers the opportunity for a direct 
comparison of the DIS results with those from e+e^ Z° exper- 
iments. The region ^ 10* GeV^, x ^ 0.50 is near the kine- 
matic limit of HERA, making it likely that data samples of about 
1000 pb^^ will be necessary for an accurate as measurement 
based on the evolution of F2. Furthermore, weak effects due 
to Z° exchange contribute to the neutral current cross sections 
for such values. It will be necessary to combine electron- 
proton and positron-proton data in order to correct for the weak 
interference terms, leading to additional possibilities for sys- 
tematic error It has been estimated [ p^ that an uncertainty on 
a^{Mz) of about 0.002 might ultimately be achieved from 
measurements of the evolution of F2 at HERA, implying a pre- 
cision of 1 .5-2.0 %. Such a precision may require a combination 
of HERA and fixed target results for F2, however [p^. 

Another possibility which has been envisioned is to oper- 
ate HERA with electron-deuteron collisions. Comparison of 
the electron-proton and electron-deuteron data would allow the 
singlet component of F2 to be extracted. A recent study [|6|] im- 
plies that this method could provide an improvement of about 
25 % in the uncertainty of as, relative to what can be achieved 
using the electron-proton data alone. 

The comments made above emphasize the relevance of 
considering electron-proton, positron-proton and electron- 
deuteron options for LEPxLHC. A detailed estimate of the as 
precision achievable using LEPxLHC has not yet been made. 
Assuming that there is not a great difference between the sys- 
tematic sources of uncertainty at HERA and LEPxLHC, it may 
be presumed that an as measurement with a precision on the 
order of 2 % is also possible at this latter facility. We note that 
LEPxLHC offers the possibihty for an accurate determination 



of as in the range of 2-3 10^ GeV^, i.e. the same range 
as a 500 GeV e+e" collider 

We therefore conclude that an as result with a precision of 
about 2 % is a possibility for HERA at a value of about 
lO"' GeV^. Extrapolating to LEPxLHC, a measurement of 
similar accuracy may be possible at a value of about 2- 
3-105 GeV^. 

Evolution ofxF^ 

The evolution of xF^ offers an advantageous method to 
measure as because it is independent of the gluon distribution 
function g{x,Q^) over the entire x range, as noted above. Mea- 
surements of xFj are best obtained using the difference between 
the 1/ and 77 cross sections for scattering on unpolarized tar- 
gets 1 15 1. These measurements require a fixed target program in 
order to collect adequate collision statistics. There is an active 
experiment at Fermilab (the NuTeV Collaboration), which is 
expected to improve the precision on a-^{l\Iz) to about 2.5 % 
within the next few years, using this method [|7[. This is likely 
to become one of the world's most precise measurements of 
as and to remain so for some time. The uncertainties on the 
NuTeV result are roughly evenly divided between statistical and 
systematic sources, with the systematic uncertainty dominated 
by imprecise knowledge of the neutrino beam flux and of the 
calorimeter energy scale. 

To improve the precision of the as result from this technique 
yet further, higher statistics from tagged neutrino beams will be 
necessary (tagged neutrino beams allow an event-by-event de- 
termination of the incident neutrino energy, as well as apriori 
knowledge of whether the interaction was caused by a neutrino 
or an antineutrino). The future facilities which could poten- 
tially provide beams for a precise xF^ measurement of as are 
therefore the primary Tevatron beam with an upgraded luminos- 
ity, such as TeV33, and the LHC. Given the good result for as 
which is anticipated from the ongoing experiment, mentioned 
above, and given the improvement in accuracy expected from 
higher statistics and the introduction of event-by-event neutrino 
tagging, it is plausible that this method can provide an as mea- 
surement with 1 % precision. A study of the precision attainable 
at a LHC fixed target experiment has not yet been performed, 
however 

In conclusion, the method based on the evolution of xF^ is 
a strong candidate to provide a 1 % measurement of a^{Mz), 
assuming that fixed target programs with tagged neutrino beams 
are available at either TeV33 or the LHC. We note that the nec- 
essary matrix elements are akeady available at NLO [^. Since 
the /3-function is known to three loops, all that is needed for a 
full NNLO analysis of as using the xF^ method are the split- 
ting functions calculated at NNLO. It is reasonable to expect 
that this result will become available and that the theoretical un- 
certainty will be below 1%. 

GLS Sum Rule 

The situation regarding the Gross, Llewellyn-Smith (GLS) 
sum rule [|l8]| is similar to that discussed above for the evo- 
lution of xFs since both methods rely on the F3 structure func- 
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tion measured in neutrino fixed target experiments. The GLS 
sum rule is based on the integral 



Jo ^ 



Table I: The estimated precision for aY^{Mz) attainable at 
future DIS experiments. 



(2) 



the QCD prediction for which has been calculated to O(as^) 
(the next-to-next-to leading order in as). This is one of the few 
quantities calculated to such a high order in QCD perturbation 
theory. The integral is evaluated experimentally using fairly 
low values of Q^, which allows small values of x. (The small 
X region is particularly important because of the 1/x weight- 
ing in (H).) The value for present experiments Jiot ] is about 
3 GeV^. Because of the low value, higher twist corrections 
are important. Furthermore, it is necessary to extrapolate into 
the unmeasured region at low x. 

Like the result for as based on the evolution of xF^, the cur- 
rent precision of the as measurement from the GLS sum rule is 
partially statistics-limited. There are several sources of exper- 
imental systematic uncertainty which are relevant for the GLS 
result and which are not relevant for the xF^ evolution result, 
however: the measurement of the absolute cross sections for 
both the and V beams, the extrapolation into the low x region, 
and higher twist corrections. The NuTeV experiment expects 
to attain a precision of about 3 % for as using the GLS sum 
rule. At TeV33 and the LHC, higher statistics, larger val- 
ues (reducing the uncertainty from higher twists) and better low 
X reach (reducing the extrapolation uncertainty) should yield 
smaller statistical and systematic errors, making a measurement 
of ay^{Mz) with a precision of 1.5 % a possibility. 

In conclusion, the GLS sum rule provides an independent 
method to determine 0^5 from a neutrino fixed target experi- 
ment, using the structure function F3. Systematic uncertainties 
should be reduced at the higher values offered by TeV33 
or the LHC, relative to the current experiments, making an as 
measurement with a precision of about 1.5 % feasible. 

Bjorken Sum Rule 

Lastly, we consider the determination of as using the Bjorken 
sum rule [n9p. The Bjorken sum rule is based on the quantity 



(3) 



where gj* and g" are the gi structure functions for proton and 
neutron targets, respectively. The Bjorken sum rule method for 
determining as differs from the others discussed here in that 
it is based on polarized cross sections. The method resembles 
the GLS sum rule technique, however, because it is based on 
an integral over a; of a structure function measured at fixed Q^, 
utilizes low measurements (for current experiments, the 
value is about 2 GeV^), requires extrapolation into the unmea- 
sured X regions, and has a QCD prediction available at the next- 
to-next-to leading order 

Li ke th e method based on the shape of F2 (item 6 in sec- 
tion [LB), it is somewhat difficult to assess the future status of 



Method 


Precision 


Facility 


evolution of F2 
at high X 


2% 


HERA, LLP X LHC 


evolution of F3 


1% 


TeV33 fixed target, 
LHC fixed target 


GLS sum rule 


1.5% 


TeV33 fixed target, 
LHC fixed target 


Bjorken sum rule 


2.5% 


Future polarized 
DIS experiments 



the as result obtainable from the Bjorken sum rule because it is 
only recently that this method has been used to determine as- 



The current result (about 4% accuracy [|12[]) is quite precise 
by current standards, however Given that additional polarized 
structure function data are currently being collected at CERN, 
DESY and SLAC, and that additional experiments are being 
planned, it can be anticipated that a reduction in the uncer- 
tainty in as from this method will be possible. Many sources 
of systematic uncertainty (higher twists, extrapolation into the 
unmeasured x region, measurement of the absolute cross sec- 
tions) are common between this method and the GLS one. The 
Bjorken sum rule measurement is complicated by its reliance on 
polarized targets and probes, however, and thus has sources of 
systematic uncertainty which are not present for the GLS mea- 
surement. Therefore, we presume that the ultimate accuracy 
for ay^{Mz) achievable from the Bjorken sum rule for exper- 
iments currently running or being planned lies between the cur- 
rent precision (4 %) and that which we estimate will be achiev- 
able from the GLS sum rule (1.5 %). Thus, an estimate of about 
2.5 % precision seems justifiable. 

Although no study has been done at this point, we wish to 
emphasize that some of the systematics which degrade the ac- 
curacy of the Bjorken sum rule measurement of as, including 
those due to higher twist and the low x extrapolation, may im- 
prove with a high statistics, high energy beam. Such a beam 
would be available if a high energy e+e" collider were con- 
structed with longitudinal polarization and a fixed-target capa- 
bility. In this way, it is plausible that the Bjorken sum rule mea- 
surement could be more accurate than the estimate given above. 

D. Conclusion for a Precise as Result from DIS 

Table || summarizes our estimates of the precision which 
might be attainable for ay^{Mz) from DIS experiments at fu- 
ture colUders. These estimates are mostly based on extrapola- 
tions from current experiments rather than on detailed studies 
of future facilities. The best prospect for a 1 % measurement 
of a-g^{Mz) from a DIS experiment is from a fixed target neu- 
trino facility at a hadron collider with high flux, tagged v and 
V beams. The most promising measurement technique is the 
observation of the evolution of xF^^. 

We again emphasize, however, the importance of accurate as 
measurements at widely different values. The DIS results 
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based on xFt, offer the possibility for a precise measurement 
of 015 in the range from about 5 to 20 GeV^. Those based 
on the high x region of F2 offer the possibiUty for an accurate 
measurement at a much larger value, up to about 10^ GeV^ 
for LEPxLHC. 



III. THE HADRON SPECTRUM 

Lattice QCD is, so far, the only systematic, first principles ap- 
proach to nonperturbative QCD. The experimentally observed 
hadron spectrum, like the high-energy observables discussed in 
the other sections, provide us with information on the free pa- 
rameter of QCD, as- Determinations of as from the experi- 
mentally observed hadron spectrum, based on lattice QCD, are 
thus complementary to determinations which are based on per- 
turbative QCD. 

While an introduction to lattice QCD is beyond the scope of 
this report (see for pedagogical introductions and reviews), 
we shall, in the following, outline the strategy for determina- 
tions of as based on lattice QCD. In general, determinations of 
as can be divided into three steps: 

The first step is always an experimental measurement. In as 
determinations based on perturbative QCD this might be a cross 
section or (ratio of) decay rates. In determinations based on lat- 
tice QCD this is usually a hadron mass or mass splitting, for ex- 
ample the mass of the p meson, or a better choice, spin-averaged 
splittings in the charmonium and bottomonium systems. In "lat- 
tice language" this step is often referred to as "setting the scale" 
(see section ^11. A ). 

The second step involves a choice of renormalization scheme. 
In perturbative QCD the standard choice is the MS scheme. 
With lattice QCD a nonperturbative scheme may be chosen, and 
there are many candidates. In order to compare with perturba- 
tive QCD, any such scheme should be accessible to perturbative 
calculations (without excessive effort). 

Finally, the third step is an assessment of the experimental 
and theoretical errors associated with the strong coupling deter- 
mination. This is of course the most important (and sometimes 
also the most controversial) step as it allows us to distinguish 
and weight different determinations. The experimental errors 
on hadron masses are negligibly small in lattice determinations 
of as at this point. The theoretical errors that are part of as 
determinations based on perturbative QCD include higher or- 
der terms in the truncated perturbative series and the associated 
dependence on the renormalization scale, and hadronization or 
other generic nonperturbative effects. In lattice QCD the theo- 
retical errors include (but are not limited to) discretization errors 
(due to the finite lattice spacing, a ^ 0), finite volume effects, 
and errors associated with the partial or total omission of sea 
quarks. 

The consideration of systematic uncertainties should guide us 
towards quantities where these uncertainties are controlled, for 
a reliable determination of as- As has been argued by Lepage 
quarkonia are among the easiest systems to study with lat- 
tice QCD, since systematic errors can be analyzed easily with 
potential models if not by brute force. 



Finite-volume errors are much easier to control for quarko- 
nia than for light hadrons, since quarkonia are smaller Lattice- 
spacing errors, on the other hand, can be larger for quarkonia 
and need to be considered. This error can be controlled by 
studying the lattice spacing dependence of physical quantities 
(in physical units). The lattice spacing is reduced (while keep- 
ing the physical volume of the lattice fixed) until the error is 
under sufficient control. The source of the lattice spacing de- 
pendence are the discretizations used in the lattice lagrangian 
(or action). Thus, an alternative to reducing the lattice spac- 
ing in order to control this systematic error is the use of better 
discretizations. This procedure is generally referred to as im- 
proving the action. For quarkonia, the size of lattice-spacing 
errors in a numerical simulation can be anticipated by calcu- 
lating expectation values of the corresponding operators using 
potential-model wave functions. They are therefore ideal sys- 
tems to test and establish improvement techniques. 

A lot of the work of phenomenological relevance is done in 
what is generally referred to as the "quenched" (and sometimes 
as the "valence") approximation. In this approximation gluons 
are not allowed to split into quark - anti-quark pairs (sea quarks). 
This introduces a systematic error into the calculation. How- 
ever, for quarkonia, a number of calculations now exist which 
partially include the effect of sea quarks, thereby significantly 
redu cing th is s ystem atic error. This is further discussed in sec- 
tions |iLA| and |lL^ 



A. Determination of the Lattice Spacing and the 
Quarkonium Spectrum 

The experimental input to the strong coupling determination 
is a mass or mass splitting, from which by comparison with 
the corresponding lattice quantity the lattice spacing, a, is de- 
termined in physical units. For this purpose, one should iden- 
tify quantities that are insensitive to lattice errors. In quarkonia, 
spin-averaged splittings are good candidates. The experimen- 
tally observed IP- IS and 2S-1S splittings depend only mildly 
on the quark mass (for masses between mi, and TOc). Figure ^ 
shows the observed mass dependence of the IP- IS splitting in 
a lattice QCD calculation. The comparison between results 
from different lattice actions illustrates that higher-order lattice- 
spacing errors for these splittings are small [23, 230. 
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Figure 2: The IP-IS splitting as a function of the IS mass (sta- 
tistical errors only) from Ref. [p3[]; □: O(a^) errors; x: 0{a) 
errors. 
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Two different formulations for fermions have been used in lat- 
tice calculations of the quarkonia spectra. In the non-relativistic 
limit the QCD action can be written as an expansion in powers 
of (or 1 / to), where v is the velocity of the heavy quark inside 
the bound state [24|; Henceforth, this approach shall be referred 
to as NRQCD. Lepage and collaborators [ p5| ] have adapted this 
formalism to the lattice regulator Several groups have per- 
formed numerical calculations of quarkonia in this approach. In 
Refs. |7|] the NRQCD action is used to calculate the bb and 
cc spectra, including terms up to 0{mv*) and 0{a^). In addi- 
tion to calculations in the quenched approximation, this group 
is also using gauge configurations that include two flavors of 
sea quarks with mass rriq ~ ^to^ to calculate the bb spectrum 



[g2| g8||. The leading order NRQCD action is used in Ref 
for a calculation of the bb spectrum in the quenched approxima- 
tion. 

The Fermilab group [ ^0| ] developed a generalization of pre- 
vious approaches, which encompasses the non-relativistic limit 
for heavy quarks as well as Wilson's relativistic action for light 
quarks. Lattice-spacing errors are analyzed for quarks with ar- 
bitrary mass. Ref. uses this approach to calculate the bb and 
cc spectra in the quenched approximation. The authors consid- 
ered the effect of reducing lattice-spacing errors from 0{a) to 
O(a^). The SCRI collaboration [ |3l| ] is also using this approach 
for a calculation of the bb spectrum using the same gauge con- 
figurations as the NRQCD collaboration with n/ ~ 2 and an 
improved fermion action (with O(a^) errors). 

All but one group use gauge configurations generated with 
the Wilson action, leaving O(a^) lattice-spacing errors in the 
results. The lattice spacings, in this case, are in the range a ~ 
0.05 — 0.2 fm. Ref. p^ ] uses an improved gauge action together 
with a non-relativistic quark action improved to the same order 
(but without spin-dependent terms) on coarse (a ~ 0.4 — 0.24 
fm) lattices. The results for the bb and cc spectra from all groups 
are summarized in Figures ^ and ^. 

The agreement between the experimentally-observed spec- 
trum and lattice QCD calculations is impressive. As indicated 
in the preceding paragraphs, the lattice artifacts are different for 
all groups. Figures |^ and Q therefore emphasize the level of 
control over systematic errors. 

Results with two flavors of degenerate sea quarks have now 
become available from a number of groups [ p2[ p3| , |4| ^8[|, with 
lattice-spacing and finite-volume errors similar to the quenched 
calculations, significantly reducing this systematic error. How- 
ever, several systematic effects associated with the inclusion of 
sea quarks still need to be studied further. They include the de- 
pendence of the quarkonium spectrum on the number of flavors 
of sea quarks, and the sea-quark action (staggered vs. Wilson). 
The inclusion of sea quarks with realistic light-quark masses is 
very difficult and can, at present, only be done by extrapolation 
from TOq ~ 0.3 — O.Stos to rriu^d- However, the dependence 
of the quarkonium splittings on the sea quark masses can be 
analyzed with chiral perturbation theory to guide the ex- 
trapolation. 
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B. Definition of a Renormalized Coupling 

Within the framework of lattice QCD the conversion from the 
bare to a renormalized coupling can, in principle, be made non- 
perturbatively. In the definition of a renormalized coupling, 
systematic uncertainties should be controllable, and at short 
distances, its (perturbative) relation to other conventional def- 
initions calculable. For example, the renormalized coupling, 
ay, can be defined from the nonperturbatively computed heavy- 
quark potential [|3^. An elegant approach has been developed 
in Ref. [^, where a renormalized coupling is defined non- 
perturbatively through the Schrodinger functional. The authors 
compute the evolution of the coupling nonperturbatively using a 
finite size scaling technique, which allows them to vary the mo- 
mentum scales by an order of magnitude. The same technique 
has also been applied to the renormalized coupling defined from 
twisted Polyakov loops |j3^]. The numerical calculations in- 
clude only gluons at the moment. However, the inclusion of 
fermions is possible. Once such simulations become available 
they should yield very accurate information on as and its evo- 
lution. A renormalized coupling can also be defined from the 
three-gluon vertex, suitably evaluated on the lattice [^9|]. 

An alternative is to define a renormalized coupling through 
short distance lattice quantities, like small Wilson loops or 
Creutz ratios which can be calculated perturbatively and by nu- 
merical simulation. For example, the coupling defined from 
the plaquette (the smallest Wilson loop on the lattice), «□ 
—3 In (Tr Un) /47r, can be expressed as 



la = ap{q)[l - (1.19 + 0.07 n/)ap(g)] 



(4) 



atq — 3.41/a, close to the ultraviolet cut-off. The coupling ap 
is chosen such that it equals ay at one-loop: 



(5) 



ap is related to the more commonly used MS coupling by 

"ms(Q) = ap(e'/'Q) (1 + f ap + C2(n/)a|, + ...). (6) 

The size of higher-order corrections associated with the above 
defined coupling constants can be tested by comparing perturba- 
tive predictions for short-distance lattice quantities with nonper- 
turbative results [^^. The comparison of the nonperturbatively 
calculated coupling of Ref. with the perturbative predic- 
tions for this coupling using Eq. (Q) is an additional consistency 
test. 

The relation of ap to a^ ,Eq. (||), has recently been calcu- 
lated to two loops [ |4l| , p3j in the quenched approximation (no 
sea quarks, Uf = 0): 



C2(n/ 0) 0.96 



(7) 



This term shifts a^^{Mz) by +0.002. Because of the unknown 
Uf dependence in the two-loop term, C2, the perturbative uncer- 
tainty is still ±0.002 (at Mz)- The extension of the two-loop 
calculation tonj 7^ will reduce this uncertainty to well below 

1% for ag^(Mz). 



C. Sea Quark Effects 

Calculations that properly include all sea-quark effects do not 
yet exist. If we want to make contact with the "real world", 
these effects have to be estimated phenomenologically or ex- 
trapolated away. 

The phenomenological correction necessary to account for 
the sea-quark effects omitted in calculations of quarkonia that 
use the quenched approximation gives rise to the dominant sys- 
tematic error in this calculation p4| |. By demanding that, 
say, the spin-averaged IP- IS splitting calculated on the lattice 
reproduce the experimentally observed one (which sets the lat- 
tice spacing, a~^, in physical units), the effective coupling of 
the quenched potential is in effect matched to the coupling of the 
effective three flavor potential at the typical momentum scale of 
the quarkonium states in question. The difference in the evo- 
lution of the zero flavor and 3,4 flavor couplings from the ef- 
fective low-energy scale to the ultraviolet cut-off, where as is 
determined, is the perturbative estimate of the correction. 

For comparison with other determinations of as, the MS cou- 
pling can be evolved to the mass scale. An average of 
Refs. [H |4| yields for as from calculations in the quenched 
approximation: 



0.110 ±0.006 . 



(8) 



The phenomenological correction described in the previous 
paragraph has been tested from first principles in Ref. [P3[]. The 
2-loop evolution of 71/ = and nj — 2 MS couplings - ex- 
tracted from calculations of the cc spectrum using the Wilson 
action in the quenched approximation and with two flavors of 
sea quarks respectively - to the low-energy scale gives consis- 
tent results. After correcting the two flavor result to = 3 in 
the same manner as before and evolving a^^ to the mass, 
they find [51 



0.111 ±0.005 



(9) 



in good agreement with the previous result in Eq. (|]). The total 
error is now dominated by the rather large statistical errors and 
the perturbative uncertainty. 

The most accurate result to date comes from the NRQCD col- 
laboration [ p2[ ^ ]. They use results for as from the bb spec- 
trum with and two flavors of sea quarks to extrapolate the 
inverse coupling to the physical three flavor case directly at the 
ultraviolet momentum, g — 3.41/a. They obtain a result consis- 
tent with the old procedure. Recently, they have begun to study 
the dependence of as on the masses of the sea quarks. Their 
preliminary result is: 



.(3)/ 



.2 GeV) = 0.195 ± 0.003 ± 0.001 ± 0.004 . (10) 



The first error is statistics, the second error an estimate of resid- 
ual cut-off effects and the third (dominant) error is due to the 
quark mass dependence. The conversion to MS (including the 
2-loop term in Eq. (^ and evolution to the Z'^ mass then gives: 



J5) 



a^(Afz) = 0.118 ±0.003, 



(11) 
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where the error now also includes the perturbative uncertainty 
from eq. (||). A similar analysis is performed in Ref. [Q| on 
the same gauge configurations but using the Wilson action for 
a calculation of the cc spectrum. The result for the coupling is 
consistent with Refs. |22, 33| ]. 

The preliminary calculation of the SCRI collaboration [ |3l[ | 
(uf — 2) can be combined with the result of Ref. [^^. Using 
the same analysis as in Ref. P4I gives 1 28 1 



by the SLD Collaboration [|46[|. Similar studies have been pub- 
lished by the LEP experiments [p7[|. The SLD result 



a!^(Mz) = 0.116 ±0.003 



(12) 



nicely consistent with Eq. (pTj). Clearly, more work is needed 

, especially in calcu- 
li. In partic- 



to confirm the results of Eqs. ( |1 1| ) and ([K 
lations with heavy quark actions based on Ref. 
ular, the systematic errors associated with the inclusion of sea 
quarks into t he simulation have to be checked, as outlined in 
section III. A 



D. Conclusions 

Phenomenological corrections are a necessary evil that enter 
most coupling constant determinations. In contrast, lattice QCD 
calculations with control over all sources of systematic error 
can, at least in principle, yield truly first-principles determina- 
tions of as from the experimentally observed hadron spectrum. 

At present, determinations of as from the experimentally 
measured quarkonia spectra using lattice QCD are compara- 
ble in reliability and accuracy to other determinations based 
on perturbative QCD from high energy experiments (see Fig- 
ure 1). The phenomenological corrections for the most impor- 
tant sources of systematic errors in lattice QCD calculations of 
quarkonia have already been replaced by first principles calcu- 
lations. This has led to a significant increase in the accuracy of 
as determinations from quarkonia. 

Still lacking for a first-principles result is the proper inclusion 
of sea quarks. A difficult problem in this context is the inclusion 
of sea quarks with physical light quark masses. At present, this 
can only be achieved by extrapolation (from niq ~ 0.3 — O.Stos 
to mu,d)- Given sufficient numerical results on the light quark 
mass dependence, chiral perturbation theory can be used for the 
extrapolation [|3^]. These calculations can most likely be per- 
formed with currently available computational resources lead- 
ing to first-principles results for the quarkonia spectra. They 
should, in turn, yield determinations of ay^{Mz) with a total 
uncertainty below 1 %. 

IV. ELECTRON-POSITRON ANNIHILATION 

The measurement of as via hadronic observables in e+e^ 
annihilation is a mature subject. Prospects for the accurate mea- 
surement of as in high-energy e+e" annihilation have been un- 
der discussion for some time [^. To assess the potential for an 
a s measurement from this method, it is useful to examine an 
experimental analysis in detail, to assess those areas in which 
the uncertainties might be reduced in the future. For this pur- 
pose, we choose a recent comprehensive study of as published 



ajjaiMz) = 0.1200 ± 0.0025 ± 0.0078 



(13) 



was derived from the consideration of 15 different infrared-safe 
hadronic observables, including various event shape parame- 
ters, jet rates derived with several different jet finding schemes, 
and energy-energy correlations. The ±0.0025 experimental 
error received contributions of ±0.0009 from event statistics, 
and ±0.0024 from detector-related uncertainties. The ±0.0078 
theoretical uncertainty resulted from contributions of ±0.0024 
from uncertainties in the hadronization process, and ±0.0074 
from missing higher orders in the perturbative calculation of the 
15 observables. Currently, all 15 observables have been calcu- 
lated to next-to-leading order in as- In addition, for six of the 1 5 
observables, the leading and sub-leading logarithms have been 
resummed and combined with the next-to-leading order calcu- 
lations. 

This breakdown of the uncertainty provides a basis for es- 
timating the accuracy of a similar measurement of as at an 
electron-positron collider of cms energy y/s — 500 GeV, such 
as the proposed NLC. Statistically, the SLD measurement was 
performed with the sample of 37,000 hadronic (e+e^ qq) 
events remaining in the 1993 SLD event sample after the ap- 
plication of hadronic event selection cuts. At a design luminos- 
ity of 5 X 10"^^ cm~^sec~^, with a Born-level cross section of 
3.1 pb, an NLC detector would collect approximately 150,000 
e+e^ — > qq events in a "Snowmass" year of 10^ seconds. The 
effects of initial state radiation and beamstrahlung, and ineffi- 
ciencies introduced by event selection (to be discussed below), 
reduce this to approximately 25,000 e+e" qq events per year 
at y/s ~ 500 GeV, adequate for a statistical precision of ±1 % 
on the value of as at that scale. A well designed NLC Detector 
calorimeter should permit a substantial reduction in the ±2 % 
detector uncertainty. 

The determination of as involves the comparison of the 
hadronic observables with parton-level perturbative calculations 
which depend upon as- The relationship between the parton- 
level calculations and the final state observables is thus obscured 
by the fragmentation process. This introduces a correction, and 
corresponding uncertainty, which must be applied to the ex- 
tracted value of as. It is generally expected that effects 
which alter the relation between the perturbative parton-level 
calculations of observables, and the actual hadron-level observ- 
ables, scale as an inverse power of the momentum transfer Q, 
so that for some observable O, 



(14) 



Typically, n equals 1 or 2. On the other hand, the perturbative 
evolution of as scales roughly as (InQ)^^. Thus, one expects 
the relative uncertainty on as due to fragmentation effects to 
scale as \nQ/Q- As a result, the ~ 10 % correction applied to 
the value of as extracted from hadronic observables at the Z*^ 
pole is expected to reduce to a ^ 2 % correction at ^/s = 500 
GeV, with an uncertainty of 1 % or less. 
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In addition to the fragmentation, the relationship between the 
perturbative parton-level calculation and the measured observ- 
ables is compromised by missing higher orders in the perturba- 
tive expansion. In the SLD analysis, this uncertainty was esti- 
mated to be l\ay^{Mz) = ±0.0074 by varying the renormal- 
ization scale of the perturbative calculation over a range per- 
mitted by consistency with the hadronic data, and observing the 
corresponding variation in the extracted value of 0:5. Current 
perturbative calculations are done to order thus, uncertain- 
ties due to missing higher orders should scale as a\, leading to 
an uncertainty of ±0.003-0.004 at = 500 GeV. Evolving 
this back to the benchmark scale Q^=M^ using the three-loop 
QCD /3-function yields an uncertainty of ±0.005-0.006, or 4- 
5 % relative, on the value of a^{Mz) extracted from hadronic 
observables at the NLC. Should next-to-next-to-leading order 
perturbative calculations become available, it should thus be 
possible to approach the target uncertainty of ±1 %. 

A final issue associated with the measurement of as in e+e~ 
annihilation at large cms energy is that of identifying a clean 
sample of e+e^ ^ QQil 7^ t) events. At ^/s = 500 GeV, with- 
out event selection cuts, qqiq 7^ production (crsorn — 3.1 pb) 
has a substantially smaller cross section than W^W" produc- 
tion {oBorn — 7.0 pb), as well as significant backgrounds from 
Z^Z^ ((TBorn = 0.4 pb) and ti (aBorn = 0.3 pb) production. A 
study performed by the European Linear Collider QCD Work- 
ing Group identified a set of kinematic cuts which select 
an 83 % pure sample of e+e^ qq (q ^ t) events. However, 
these cuts substantially impacted the hadronic distributions of 
the remaining qq events, leading to large 20 %) corrections 
to the extracted value of as- To this end, in preparation for 
the Snowmass Workshop, a Monte Carlo study [ 5^ was under- 
taken in which one of the European Working Group kinematic 
cuts - the requirement that at least one hemisphere have a recon- 
structed invariant mass of less than 13 % of the visible energy 
in the event - was removed. Instead, events were used only if 
they were produced with the right-handed electron beam (to re- 
move W^W^ background), and if they gave no indication of 
the presence of B hadrons in the vertex detector (to eliminate tt 
background). For an electron beam polarization of Pg = 80 % 
(Pe = 90 %), this yielded an 82 % (87 %) pure "Snowmass" 
sample of e+e^ qq (q ^ t) events. A comparison of 3-jet 
rates between a pure Monte Carlo sample of e+e^ ^qqiq 
events, and the sample identified by the Snowmass cuts, indi- 
cated that corrections due to the Snowmass event selection are 
substantially less than 5 %. Thus, with these cuts, the uncer- 
tainty on as due to the event selection process is expected to be 
well within the target of ±1 %. It should be noted that electron 
beams with 80 % polarization, and bunch populations exceed- 
ing that required for the operation of the NLC, are already in 
use at the SLAC Linear Collider, and that polarized running is 
part of the base-line proposal for the NLC [pl|. 



As a final note, it has been pointed out [ p2| ] that the high lu- 
minosity of an e+e^ linear collider, combined with the rise in 
the e+e^ qq cross section with falling y/s, may make it feasi- 
ble to precisely constrain the evolution of as over a wide range 
of in a single experiment. The execution of such a pro- 



gram would have an impact on the design of the high energy 
e+e~collider. 

Thus, e+e~ annihilation at high energy appears to be a 
promising avenue towards the measurement of a^{Mz) to a 
relative uncertainty of ±1 %. Furthermore, the high momen- 
tum transfer scale associated with the measurement (Q^ ~ s — 
(500 GeV)^) makes this approach an important component of 
the program to constrain the possible anomalous running of as- 
For this accuracy to be achievable, next-to-next to leading order 
(0{ag)) calculations of e+e^ event shape observables will be 
required. 

V. pp (pp) COLLISIONS 

The greatest potential to extend measurements of as to large 
values of the momentum transfer scale resides with the pp 
(pp) colliders. In adition, many approaches to the measurement 
of asin pp (pp) collisions produce a range of values for as over 
a broad lever-arm in Q^- For example, the inclusive jet Et spec- 
trum from the Tevatron extends out to almost 500 GeV (see Fig- 
ure ^ [^^), providing sensitivity to 0^5 over a range in momen- 
tum transfer extending from 50 GeV to values nearly equivalent 
to that proposed for the next generation of electron-positron and 
electron-proton colliders. The LHC, currently scheduled to be- 
gin delivering beams in the middle of the next decade, will ex- 
tend this reach to several TeV 
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Figure 5: The preliminary CDF and DO Run lb data compared to 
NLO QCD using CTEQ4M parton distributions. Experimental 
points normalized as indicated. This figure is reproduced from 
Reference 



On the other hand, measurements of 0^5 from hadron collid- 
ers have not yet approached the level of accuracy achieved by 
the most accurate approaches. For example, a typical approach 
to constraining 0^5 in pp collisions is to study the ratio of 14^ -1- 1 
jet events to -1- jet events [|54|]. Experimental systematics, 
such as energy scale and resolution uncertainty, introduce large 
errors (± ^ 0.015) in the value of as extracted from this ra- 
tio. In addition, since gluons liberated from the nucleon sea can 
themselves form jets, the measurement is sensitive to the parton 
distributions used in calculating the H- 1 jet rate. For the most 



10 



recent measurement of this ratio [p5[|, the DO collaboration did 
not report a value for as, because of an inconsistency between 
the best fit value of as used in the parton distribution function 
and that used in the perturbative matrix element. 

Current work on the measurement of as with pp data thus 
concentrates on developing approaches which remove or reduce 
the sensitivity of the method to variations in the parton distribu- 
tion functions, and to experimental parameters. This work is 
still in its early stages, but a number of promising ideas are be- 
ing pursued. 

For example, a generalization of the W + jet method is the 
measurement of the jet cross section ratios 



R„ 



a{V +{n+ l)jets) 
a{V + njets) 



(15) 



where V = W^, is a vector boson. The UA2 and DO mea- 
surements involved Rq; for n ^ Q, however, the contribution 
from sea gluons, and thus the dependence on the parton distri- 
bution functions, largely cancels in the ratio. Another approach 
that is being pursued is the measurement of the p± spectrum of 
production - for this measurement, there is no dependence 
on experimental errors such as hadronic energy scales and reso- 
lution, jet algorithm definitions, or hadronization, although the 
measurement still requires apriori knowledge of parton distri- 
bution functions, and only measures as at the single momen- 
tum transfer scale Q^=M^. Finally, fits to the triple differential 
cross section pq] 
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are being explored, which can simultaneously constrain as, the 
gluon distribution function fg (x), and the non-singlet quark dis- 
tribution function f2{x), thus removing the uncertainty due to 
poorly constrained parton distribution functions. 

All of these studies have only recently been started [|7[, in- 
spired by the large data samples available with the completion 
of Tevatron Run 1. Thus, it will be several years before the po- 
tential for the measurement of 0:5 at hadron colliders is fully 
understood. Finally, it should be noted that this method (like 
any other) requires that at least NNLO perturbative calculations 
be completed for a determination of a^{Mz) with 1 % accu- 
racy. However, most of the matrix elements needed here are 
identical to, or are limiting cases of, matrix elements necessary 
for NNLO calculations of hadronic observables in e+e~ anni- 
hilation. 

VI. CONCLUSIONS 

Table || lists the methods we consider promising for accurate 
as determinations. The items listed above the double line are 



established methods for as measurements, and we can evalu- 
ate their potential for 1 % accuracy with reasonable confidence. 
The items listed below the double line are either expected to 
yield somewhat less accurate determinations of as, or they are 
less established methods which need further study to better eval- 
uate their potential for as determinations with 1 % accuracy. 

Traditionally, DIS at relatively low momentum transfer has 
produced precise determinations of as. In particular, measure- 
ments of as via the evolution of xF^ and the GLS sum 
rule are expected to each achieve experimental accuracies of 
2-3 % in the upcoming run of the NuTeV Experiment, limited 
primarily by sample statistics, the uncertainty in the calorimeter 
energy scale, and the understanding of the composition of the 
incident neutrino beam. Thus a high flux tagged neutrino beam- 
line derived from the full energy Tevatron primary beam, and 
eventually one of the LHC primary proton beams, is a strong 
candidate for a facility which will produce a 1 % measurement 
of as at low Q^. At present, however, such a facility is not 
part of the future program of either laboratory. We also wish 
to mention the approach of measuring as via the Bjorken sum 
rule in polarized deep inelastic scattering. It is a relatively new 
method, but could yield a result as accurate as 2-3 %. Certain 
systematic limitations, such as the corrections for higher twist 
and the extrapolation of gi into the unmeasured region at low 
X, may be less problematic if polarized high energy NLC elec- 
tron beams are available for fixed target physics. This issue is 
worthy of further study. 

Lattice QCD calculations have matured considerably in the 
last few years. First principles calculations of the simplest 
hadronic systems, like quarkonia spectra, should be possible 
with current technology and computational resources. This 
implies the potential for very accurate determinations of as 
at relatively low from experimental measurements of the 
hadron spectrum, with systematic uncertainties largely indepen- 
dent of all other approaches discussed here. Since the present 
experimental errors contribute much less than 1 %, no future 
experimental facilities are required for a 1 % determination of 
ay^{Mz) from the hadron spectrum. 

e+e~ annihilation experiments at high center-of-mass ener- 
gies are promising for accurate determinations of as from mea- 
surements of jet rates and other jet variables. Such experiments 
could be performed at an NLC collider with Ecm = 500 GeV 
or higher. 

ep scattering experiments at HERA which measure the struc- 
ture function F2 at high x over a wide range of can poten- 
tially yield determinations of a-^{Mz) with about 2 % accu- 
racy. It should be noted that a future LEPxLHC facility can 
potentially probe momentum transfers of ^ (500 GeV)^. 

Experiments at hadron colliders (pp or pp) have the high- 
est potential energy reach. Accurate determinations of as will 
require either concurrent extractions of the parton distribution 
functions (PDFs) from the same experiment or prior knowledge 
of the PDFs (with error bars) over the range of x probed by the 
process under study. Feasibility studies are underway now for 
the Tevatron experiments, and it is expected that the potential 
for providing an accurate measurement of as in high energy 
hadron collisions will be understood within the next few years. 
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Table 11: Summary of methods for potential 1 % determinations of aj^{Mz)- The methods listed below the double Une are either 
considered to yield somewhat less accurate determinations (F2 at high x at HERA, Bjorken sum rule), or they have not yet been 
fully established (pp and pp collisions). 



Process 


Approach 


NNLO Calculation 


Energy Scale 


FGcilities 


DIS 


evolution of F3 


Partial 


2-15 GeV 

Z-^-j \Jc V 


TeV33 fixed target 
LitiK^ nxea target 


DIS 


GLS sum rule 


Available 


few GeV 
^ i u vje V 


TeV33 fixed target 
L n v-^ uxeu Large I 


Hadron spectrum 


Spin-averaged hh and cc sphttings 


lattice QCD 


few-lOGeV 


none 




Hadronic observables 


Partial 


500 GeV 


NLC 


Polarized DIS 


Bjorken sum rule 


Available 


few GeV 
- 10 GeV 


SLAC, DESY, HERA 
NLC fixed target 


DIS 


evolution of F2 at high x 


Partial 


few-100 GeV 
< 500 GeV 


HERA 
LEPxLHC 


pp 
pp 


Jet properties 


Partial 


100-500 GeV 
< few TeV 


Tevatron 
LHC 



The reach of the LHC is substantially larger than that of 
any other accelerator mentioned here. Should it prove possible 
to accurately determine as in hadronic collisions, the construc- 
tion of a higher energy collider would extend this reach even 
further. 

In summary, the goal of measuring as to an accuracy of 1 %, 
with a number of complementary approaches, and over a wide 
range of Q^, seems feasible. The complete program will likely 
require a number of new facilities. At low Q^, to approach a 
precision of 1 % in DIS experiments, it will most likely be nec- 
essary to establish a tagged neutrino beam facility utilizing the 
full energy Tevatron beam, or eventually one of the LHC pro- 
ton beams. The determination of as from the hadron spectrum 
using lattice QCD is the only method without facility implica- 
tions. At high Q^, a measurement of as in e+e" annihilation 
fits quite naturally into the physics program for the proposed 
Next Linear Collider The potential for complementary as de- 
terminations in pp collisions at the Tevatron (and pp collisions 
at the LHC) still needs further study, as do measurements of F2 
at high X at HERA or at a possible LEPxLHC lepton-hadron 
collider. 
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